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Controlled release systems composed of hydrophobic polymers such as ethylene-vinyl 
acetate copolymer have proven useful for releasing various polypeptides and other macro- 
molecules for over 100 days. However, the release mechanism has never been elucidated 
evidence by microscopy is presented suggesting that release occurs through intercon- 
necting pores formed by the macromolecules themselves. A mathematical model has been 
developed and used to predict the release rates of different proteins 



INTRODUCTION 

Biocompatible polymers, such as ethylene- 
vinyl acetate copolymer (EVAc), permit the 
controlled release of macromolecules [1]. 
These polymer systems have a wide variety 
of applications. For example, EVAc polymers 
have been used as release systems for chemo- 
tactic [2,3] and growth factors [4,5], as 
components of bioassays for informational 
macromolecules [6,7] and histochemical 
markers in neurological studies [8,9], and 
also as delivery systems for insulin [10], 
interferon [11], and antigens [12]. How' 
ever, the mechanism by which macromolec- 
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NY, U.S.A. 

**Present address: School of Pharmacy, University of 
California, San Francisco, CA 94143, U.S.A. 
***To whom correspondence should be addressed at 
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ular release occurs has not been elucidated. 
The polymers used are impermeable to 
molecules larger than 600 daltons; nonethe- 
less, complete release of substances as large 
as 2 X 10 6 daltons for over 100 days has been 
demonstrated [1], We now report that the 
incorporation of macromolecules into the 
normally non-porous polymer matrix results 
in formation of a tortuous, interconnected 
pore network. It is suggested that diffusion 
of the macromolecules through this network 
provides the basis for controlled release. 



MATERIALS AND METHODS 
Kinetic studies 

Bovine serum albumin (BSA) (M w = 69000), 
/J-lactoglobulin A(M W = 37000), and lysozyme 



0l68-3659/85/$03.30 © 1985 Elsevier Science Publishers B.V. 



-127 



169, 



177 



183 



se 



191- 



1971 



■4 

205: 



217 



225 I 



2331 



.239. 



.247 



.2491 



uced rates are. 
. India. Israej^ 
r copies by a 

i Netherlands 

t be honoured) 

sevier Science.! 
12) 916-12501. 



journal of 



controlle 
releas 



OFFICIAL JOURNAL OF THE CONTROLLED RELEASE SOCI 





I 



Elsevier 



260 

2* 

(M w = 14000) (all from Sigma Chemical Co., 
St. Louis, MO) ppwders were sieved into 
specified particle size ranges [13]. Each pro- 
tein powder was dispersed in a 10% (w/v) 
ethylene— vinyl acetate (Elvax 40P, DuPont, 
40% by weight vinyl acetate, washed with 
solvents to remove impurities [14]) solution 
in methylene chloride, and the suspensions 
were cast in a flat glass mold at -80°C. The 
low temperature caused the mixtures to 
congeal quickly, preventing drug migration. 
The resulting sheets were dried in two 48- 
hour stages, first at -20°C, and then under 
vacuum at 20°C [13]. 

To test release kinetics, nine 1 cm X 1 cm 
X 1 mm slabs were cut from each sheet, and 
coated on five faces using paraffin [13], 
leaving one 1 cm X 1 cm face exposed. 
Two straightened stainless steel autoclips 
were then pressed into the paraffin on the 
back of each slab to anchor it down when 
placed in the release medium. The release 
medium (0.9% NaCl solution) was placed 
in 10 ml amounts into 20 ml scintillation 
vials. The slabs were placed on a shaker 
as described previously [15]. At each time 
point, slabs were moved to vials containing 
fresh saline, and the old solutions were 
spectrophotometrically (280 or 220 nm) 
assayed for protein content [13] . 

Protein particle densities were determined 
using a pycnometer with methylene chloride 
as the solvent. Before release, the porosity 
of a slab was determined by dividing the 
protein concentration in the slab by the pro- 
tein particle density. At the end of the re- 
lease experiment, porosity was again assessed 
by liquid leaching of salicylate [15] . In 
general, the porosity values before and 
after release agreed to within 5%. Thickness 
of the polymer slabs was measured using a 
micrometer (Stamet Co., Athol, MA). The 
standard deviation for thickness measure- 
ments (n = 8) was less than 4% for each slab. 

Microscopic studies 

For optical microscopy, the following 
procedure was used: 10 nm thin sections 



were obtained using a cryogenic microtome 
set at -25°G, together with a microtome 
knife (Cryo-Cut Cryostat Microtome, Model 
845 with a 4%" microtome knife Type 
942, American Optical Corp., Scientific 
Instruments Division, Buffalo, NY). One 
milliliter of embedding medium (Tissue-Tek 
II O.C.T. Compound 4583, Lab-Tek Products 
Division, Miles Inc., Naperville, IL) was 
poured onto the chucks. After 30 seconds, 
it froze into an opaque solid, which was 
then planed by cutting thin sections off 
its top surface. Next, about 3 mm X 1 mm 
X 1 mm of release matrix was excised from 
various locations in the original 1 cm X 1 cm 
X 1 mm slab. This piece of matrix was then 
placed on the planed embedding medium such 
that one of the 3 mm X 1 mm faces cor- 
responding to a cross-section of the matrix 
was in contact with the planed surface. 
The matrix so oriented was then buttressed 
by more embedding medium which also 
hardened after 20 seconds in the machine. 

Sections of 10 iim in thickness were then 
cut, and stuck of their own accord to the 
knife. These were retrieved by contact with 
a glass slide at room temperature. The final 
sections thus obtained were 3 mm X 1 mm 
X 10 jum, the 1 mm representing the original 
depth of the matrix. 

For observation under scanning electron 
microscopy (SEM), the following procedure 
was followed: The slabs, having been ex- 
posed to the release Medium (0.9% saline), 
were dried in order to enable the high vacuum 
conditions necessary for SEM. To prevent 
deformation during drying of the pliable 
wet slab matrices, a critical point drying 
machine (Model 11-120A, Balzers Union, 
Liechtenstein) was used. The wet slab was 
placed in the chamber of the machine and the 
water in the sample replaced by 100%ethanol. 
The ethanol was then replaced by 100% 
amyl acetate, which is miscible with liquid 
carbon dioxide. The chamber was then 
cooled to 4°C and filled with liquid carbon 
dioxide. Carbon dioxide vapor was slowly 
exhausted to air while refilling the chamber 
with more carbon dioxide liquid to remove 
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drying chamber was then heated to 40°C 
and the pressure increased to 85 atm As 
soon as this temperature and pressure were 

r q uTd ve to (15 ~ 20 r } ' the cha ^ e 

liquid to vapor phase of carbon dioxide 
occurred, and the sample was dry The 
pressure and temperature were then de- 
creased to ambient conditions*. The dried 
samples were then coated using an Economy 
Coater Type CVE-15 (Consolidated Vacuum 
Corporation, Rochester, NY) under vacuum 
with carbon and gold. The dried and coated 
samples were stored in tightly capped vials 
which were placed in a desiccator. 

Drug distribution studies 
t 

Slabs containing BSA (50% loading) were 
released into 100 ml volumes of 0.9% NaCl 
At various timepoints, a slab was removed 
and frozen on dry ice to terminate release 
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A eryomicrotome (Damon/IEC, -25°C) was 
used to serially cut 10 jum sections from the 
partially released polymer matrices. Four 
samples, representing four consecutive 10 
w f re th *n released into 10 ml 
0.9% NaCl solution for three days The 
release medium was then filtered to remove 
wax polymer, and embedding medium 
Protein concentration was determined by 
UV spectroscopy at 220 nm [13] . 

RESULTS AND DISCUSSION 

Viewed under an optical microscope, 
EVAc films cast without proteins appear 
as non-porous sheets (Fig. l a ). Slabs cast 
with proteins and sectioned prior to release 



fo?™!.™ P r ° CedUreS f ° r P re P a "ng EVAc slabs 
for SEM have been developed recently [16]. 
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display areas of either polymer or protein 
(Fig. lb). Slabs initially cast with proteins 
and released to exhaustion (5 months) appear 
as porous films (Fig. Id). Pores with diam- 
eters as large as 100 nm, the size of the 
sieved particles, were observed. The structures 
visualized were also confirmed by Nomarski 
(differential interference contrast) micro- 
scopy. It appears that although pure EVAc 
is impermeable to macromolecules [1], 
molecules incorporated in the matrix dissolve 
once water penetrates the matrix and are 
then able to diffuse to the surface through 
pores created as the particles dissolve. 

The scanning electron micrographs in 
Fig. 2 show that the pores are intercon- 
nected through narrow passageways. The 
passageways are necessary to permit move- 
ment of the macromolecules between pores* 

Changes in pore structure over time were 
investigated. Sections were prepared from 
matrices in the process of release (Figs, 
lb— d). We observed that (1) the pore struc- 
ture changes minimally as a function of time, 
(2) after 16 or 40 hours there is no evidence 
of a receding interface between dissolved and 



dispersed drug (Figs. lc,d), and (3) none of 
the drug remains undissolved at 40 h (30% 
release). Observations (2) and (3) differ 
from those reported for less soluble low 
molecular weight drugs such as certain ster- 
oids [17], and are probably due to the high 
solubilities of many proteins such as BSA 
(solubility >500 mg/ml) [18] . 

Figure 3 shows kinetic curves for BSA 
slabs of various drug particle sizes and load- 
ings. The captions indicate the corresponding 
measured porosities. 

We have made a number of assumptions 
in the development of a model: (1) The 
rate-limiting step for transport is drug dif- 
fusion through pores (other steps such as 
water penetration into the pores and drug 
dissolution occur in less than 40 h (Figs, lc, 
d)). (2) The effect of concentration depen- 
dence on the drug diffusion coefficient 
is not significant [19]. (3) No drug dif- 
fusion occurs through the polymer backbone 
[1]. (4) The pores are interconnected (Figs, 
lb— d, 2a,b), the porosity is uniform, and 
pore size changes minimally with time (Figs, 
lb— d). (5) The initial drug distribution is 




5P.Outer Pore 




Pore Body 



^Connecting Channel 



Fig. 2. Scanning electron micrographs (SEM) of controlled release polymers. 

(a) Surface of drug-containing EVAc matrix after termination of release (60 h). The average pore size is 101 ± 33 
/im, in the same size range as the pores and particles in Figs, lb— d. Controls consisting of pure EVAc matrices 
and EVAc matrices containing drugs before release show no such pore structure [25]. Loading = 0.25, particle 
size = 63— 149 ^m. 

(b) 10 X magnification of one of the outer pores in (a). Notice that there is a channel leading to an inner pore at 
the base of the outer pore. The pore body is the inside of the outer pore. 
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uniform (Figs. 1,6). (6) No boundary layer 
effects exist*. (7) Infinite sink conditions 
exist**. (8) Minimal effects exist due to 
osmosis or charge interaction of the drug 
with the polymer***. (9) Release occurs 
from only one face of the slab, since the 
other five faces are coated with an imper- 
meable waxt. 

With these assumptions, release is modeled 
by Fick's second law [20] : 



be _ 3 2 c 

dl~ De Jx 7 



0< x < L 



with the boundary conditions 

c(x=L,t) = 0 

and 

dc 



dx 



= 0 



x = 0 



(1) 



(2) 



(3) 



* This was verified by stirring, which would have 
disrupted boundary layers had they been present. 
Release rates of slabs stirred in containers at 2000 
rpm were compared to those that were on the shaker 
and those that were not shaken at all. Over a 400- 
hour time period (through 60% release; data were 
taken at 17 different timepoints) there was no dif- 
ference in any of the release rates. This indicates the 
lack of boundary layer effects. 

**The volume of the release medium was approxi- 
mately 100 times the volume of the slab. Increasing 
the release medium volume does not alter measured 
release kinetics. 

***Consoriant with this assumption, we found no 
effect on release rate due to increasing the ionic 
strength of the medium from 0 to 1 M NaCl [11]. 
"•"The impermeability of the paraffin wax has been 
verified by control experiments where all faces 
are coated with wax. 

^In this paper concentrations are expressed in terms 
of the whole matrix including both the aqueous 
pore space and the polymer. This is in contrast to 
the more common usage, in which concentration 
is in terms of the aqueous pore space alone. This 
allows us to write FickV second law (eqn. 1) without 
reference to the matrix porosity. 

have chosen cm 2 /h instead of cm 2 /s as the unit 
for the diffusion coefficients because it reflects the 
time scales over which observations were made. 



(corresponding to the situation where drug 
flow into the release medium occurs at 
x = L), and the initial condition 



c(xJ = 0) = C o , 0<x<L 



(4) 



In eqns. (1)— (4), L is the thickness (cm) 
of the slab, t is time (h), x is a position 
within the slab (cm), c(xj) is the local drug 
concentration (mg/cm 3 slab)* 1 " 1 ", and C 0 is the 
initial drug concentration. D e is the effective 
diffusion coefficient (cm 2 /h)^, which is 
defined and discussed below. 

The solution to eqns. (1)— (4) is [20] 



4C 0 ~ ("IT 
c(xj) = £ 

7T n = 0 2n + l 



(5) 



X exp [- (2n+l ) 2 tt 2 D e t/ 4L 2 ] cos 



(2n + l )irx 
2L 



The cumulative fraction of drug released (i.e., 
the amount released divided by the amount 
originally incorporated), is given by [20] 



M f 8 00 



= 1- 

Moo 7T 



n = 0 (2n + iy 



Xexp [-(2n+l) 2 7r 2 Z> e t/4L 2 ] 



(6) 



Figure 3 contains fits of eqn. (6) to release 
kinetic curves for BSA released from EVAc 
slabs, with various drug loadings and particle 
sizes. The effective diffusivity, D e , is a free 
parameter whose value for each loading 
and particle size is determined by the fit to 
the corresponding kinetic curve. 

We assume that D e is composed of two 
factors. The first factor is D 0 , the bulk dif- 
fusivity of. the drug molecule. The second 
factor, which we shall call F, accounts for 
geometric effects of the pore structure, 
which include tortuosity, dead-end pores, 
and constrictions between pores. 

Thus 



D e =D 0 F 



(7) 



Since D 0 is either measureable or obtainable 
from the literature, the pore structure depen- 
dent factor F can be extracted from the fits 
in Fig. 3. 
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A log— log plot of F versus porosity (see 
Fig. 4) was well fit by the function 



log 10 F = 0.463 + 5.64 log 10 e 



(8) 



where e is the porosity. Knowing this equa- 
tion for F, and using eqn. (7), we can then 
write 



D e =£> 0 (2.904e 5M ) 



(9) 



and this value of D e can be substituted into 
eqn. (6) to predict release kinetics for slabs 
with other porosities. 

A test of our model is to cast slabs using 
other proteins, measure the parameters L, 
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e, and Z) 0 , and see whether the release kinetics 
follow eqn. (6) with D e given by eqn. (9). 
This has been done for 0-lactoglobulin and 
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Fig. 4. Log-log plot of factor F = D e /D 0 as a func- 
tion of porosity for BSA matrices whose kinetics 
are shown in Fig. 3a— c. D Q for BSA is 2.52 X 10" 3 
cm 2 /h [18] (value corrected for T = 25°C). For 
each combination of loading and particle size, D was 
determined using the best fit of eqn. (6) to theWa. 
Line is best fit of data, and determines eqn. (8). 



lysozyme (Fig. 5). The solid lines are predic- 
tions based on eqns. (6) and (9), which show 
general agreement with the data. The dif- 
ferences that are observed between prediction 
and experiment may be due to contributions 
due to differences in the shapes of the protein 
powders incorporated into the device (which 
could affect pore geometry), or to deviations 
from assumptions 1—8 for the cases studied. 

An additional test of the model is to 
determine whether it can predict the time- 
dependent concentration profile of the drug 
within the matrix. Concentration profiles of 
a BSA slab with loadings, measured at four 
timepoints, are shown in Fig. 6. The curves 
in Fig. 6 are the predictions of eqn. (5), 
with D e determined by fitting eqn. (6) to 
the accompanying release kinetic data. 
The fits are quite good, considering the 
difficulties in the experimental techniques 
used to obtain the concentration profile. 

The above equations have several limita- 
tions. It is not expected that they will apply 
(a) at low loadings (<15%), where the pore 
space may not be interconnected [21], 



Beta- Lactoglobutin 75-150 pm 




Lysozyme 106-150 
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Sqrt time (h 1 ^) 
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Fig. 5. Release kinetics of /Mactoglobulin and lysozyme. Lines are predictions based on eqns. (6) and (9) Each 
data point represents the mean of 8 values. 

(a) Mactoglobulin: particle size = 75-150 »m y D 0 = 2.82 X 1(T 3 cm 2 /h [26] (value corrected for T = 25°C) 
Mandard deviations at all points are <0.04: - - loading = 0.25, e = 0.21; * - loading = 0.40 e = 0 27- ♦ - 
loading = 0.50, e = 0.32. ' 

(b) Lysozyme: particle size = 106-150 *im, D 0 = 3.74 X 10~ 3 cm 2 /h [27] (value corrected for T= 25°C). Stan- 
dard deviations at all points are <0.04: a — loading = 0.40, e = 0.27. 
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Concenirotion profile 

600r 




Fig. 6. Drug distribution inside the matrix as a 
function of time and normalized (x/L) position. 
Slabs containing BSA at loading 0.5 were released 
into 100 ml volumes of 0.9% NaCl for varying periods 
of time. A cryomicrotome (Damon/IEC) was used 
to serially section the partially released polymer 
matrices. Protein concentration was determined by 
UV spectroscopy. Data are indicated by symbols. 
Each symbol represents the mean concentration 
of 4 samples. Lines are predictions of eqn. (5). D e 
was determined directly by fitting eqn. (6) to the 
release data. C 0 was directly determined to be 484 
mg/(cm 3 matrix) using weight of protein fraction 
and volume of matrix. ■ — 2 hours; * — 6.5 hours; 
♦ — 29.5 hours; • — 172 hours. 

or (b) for macromolecules whose solubility 
is less than 250 mg/ml, where the assumption 
of rapid drug dissolution may be incorrect. 
In the latter case, moving-zone models may 
be more applicable [22] . 

It should also be recognized that the rela- 
tion between the F factor and the porosity 
(eqn. 9) is empirical, and may not be useful 
outside its known range of validity (0.15 < e 
< 0.35). 

The diffusion equations used in the current 
study are simplications of more complex 
processes. The factor F must take into ac- 
count those matrix pore geometric factors 
contributing to decreases in diffusion rates. 
Such factors may include pore " tortuosity", 
dead-end pores, and constrictions between 
pores [21]. The understanding of such 



factors will be important in further develop- 
ing macromolecular delivery systems. 

Although the present study reports in vitro 
release data, we have previously shown that 
in vitro and in vivo release rates of macro- 
molecules from identical EVAc slabs are 
identical [23]. While the release rates from 
the slabs decrease with time, constant release 
is attainable from EVAc systems with ap- 
propriate geometric design [24] . 

The present study helps to explain why 
macromolecules can slowly permeate through 
normally impermeable polymers. The data 
should be useful in the design of release 
vehicles for various polypeptides, polysac- 
charides, and other bioactive agents now 
produced by genetic engineering. Such 
substances often possess very short in vivo 
life times (e.g., growth hormone, interferon) 
[21], and conventional dosage forms cannot, 
in general, be used to deliver these drugs. 
The methodology developed here may be of 
value in the design of systems for long-term 
delivery of macromolecules. 
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